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Abstract 

The density of the Fisher zeroes, or zeroes of the partition function in the complex 
temperature plane, is determined for the Ising model in zero field as well as in a pure 
imaginary field i7r/2. Results are given for the simple-quartic, triangular, honeycomb, 
and the kagome lattices. It is found that the density diverges logarithmically at points 
along its loci. 

1 Introduction 

In the analyses of lattice models in statistical mechanics such as the Ising model, the parti- 
tion function is often expressed in the form of a polynomial in variables such as the external 
magnetic field and/or the temperature. Since properties of a polynomial are completely 
determined by its roots, a knowledge of the zeroes of the partition function yields all ther- 
modynamic properties of the system. Particularly, if the zeroes lie on a certain locus, a 
knowledge of its density distribution along the locus is equivalent to the obtaining of the 
exact solution of the problem. 

For the Ising model with ferromagnetic interactions, we have the remarkable Yang-Lee 
circle theorem |]1| which states that all partition function zeroes lie on the unit circle |2;| = 1 
in the complex z = e^^ plane, where L is the reduced external magnetic field [wesetkT = 1). 
However, the density of the Yang-Lee zeroes on the unit circle, a knowledge of which is 
equivalent to solving the Ising model in a nonzero magnetic field, is known only for the Ising 
model in one dimension. 

Fisher p| has proposed that it is also meaningful to consider partition function zeroes in 
the complex temperature plane. Indeed, he showed that for the zero-field Ising model on the 
simple quartic lattice with nearest-neighbor reduced interactions K, the partition function 
zeroes lie on two circles 

|tanhir± 1| = (1) 

in the thermodynamic limit. He further showed that the known logarithmic singularity of 
the specific heat follows from the fact that the density vanishes linearly near the real axis. 
Subsequently, the Fisher loci has been determined for the infinite triangular lattice 0], and 
for finite simple-quartic lattices which are self-dual [Q. Stephenson has also evaluated the 
density distribution on the circles in terms of a Jacobian. However, the explicit expressions 
of the density function of the Fisher zeroes do not appear to have been heretofor evaluated. 

In this paper we complete the picture by evaluating the density function. We deduce 
the explicit expressions for the density of Fisher zeroes for the simple-quartic, triangular, 
honeycomb, and kagome lattices. Density of the Fisher zeroes for the Ising model in a pure 
imaginary field L = i7r/2 are also obtained. 
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2 The simple-quartic lattice 



It is well-known that the bulk solution of spin models with short-range interactions is in- 
dependent of the boundary conditions. For the Ising model on the simple-quartic lattice, 
we shall take a particular boundary condition introduced by Brascamp and Kunz |^ for 
which the location of the Fisher zeroes is known for any finite lattice. This permits us to 
take a a well-defined and unique bulk limit, thus avoiding a difficulty encountered in the 
consideration of the Ising model on a torus 

Consider an M x 2N simple-quartic lattice with cylindrical boundary conditions in the 
direction and fixed boundary conditions along the two edges of the cylinder. The 2N 
boundary spins on each of the two edges of the cylinder have fixed fields ■■■ + + + + + + ■■■ 

and ■■■H 1 1 — respectively. This is the Brascamp-Kunz boundary condition 

Brascamp and Kunz showed that the partition function of this Ising model is precisely 



h,MK)=2''''' n n 

l<i<N l<j<M 



1 + z — z(cos 6i + cos ( 



where 

2 = sinh2i^, ei = {2i-l)Tx/2N, (pj = jtx / {M + 1). (3) 
The per-site "free energy" in the bulk limit is then evaluated as 

f = lim In Zm2n(K) 



- ln(42) + ^ / dO f d(f) ln[z + z - (cos 6 + cos 0)] 

Z OTT J —TV J —TV 

1 1 ("^ 

-ln(42;)H / dO rf0 ln[z + - 2 cosm cost;], (4) 

2 27r^ JO JO 



where u = {6 + (j))/2,v = {6 — (j))/2 and we have made use of the fact that the integrands 
are 27r-periodic. 

The partition function @ has zeroes at the 2MN solutions of 

z + z''^ = cos 9i + COS (pj, I <i < N;l < j < M. (5) 

The following lemma and corollaries are now used to determine the loci of the zeroes: 

Lemma: The regime —2 < z + z~^ < 2 of the complex z plane, where z + z~^ = real, is the 
unit circle l^l = 1. 

Proof: The Lemma follows from the fact that, by writing z = re*^, we have 

z + z~^ = H — ^ cos 9 + i(^r ^ sin 9, (6) 

so that z + z~^ = real implies either r = 1 or 6' = integer xvr. In the latter case we have 
\z + z^'^l = |r + r^^\ > 2, which contradicts the assumption, unless r = 1. It follows that we 
have always r = 1, or \z\ = 1. Q.E.D. 

Corollary 1: The regime —a < z + z^^ < b, where a,b > 2, z + z~^ = real, of the complex 
z plane is the union of the unit circle \z\ = 1 and segments z_{—a) < x < z^{—a) and 
z_{b) < X < z+{b) of the real axis, where z±{b) = (6 ± \/b'^ — 4)/2. 
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Corollary 2: The regime —a < z + z ^ < b, where a,b > 2, z + z ^ = real, of the complex z 
plane, is the regime |if | = 1 in the complex w plane, where w is the solution of the equation 

w + vj-^ = ^-(z + z'^ + %-^). (7) 
a + b\ 2 J ^ ^ 

Corollary 1 is established along the same line as in the proof of the lemma, and Corollary 
2 is a consequence of the lemma since, by construction, we have —2 < w + < 2. 

Returning to the partition function (^), since the right-hand side of is real and lies 
in [—2,2], it follows from the Lemma that the 2MN zeroes of @ all lie on the unit circle 
|sinh2i^r| = 1, a result which can also be obtained by simply setting the argument of the 
logarithm in the bulk free energy equal to zero. The usefulness of this simplified procedure 
has been pointed out by Stephenson and Couzens ^ for the Ising model on a torus. But 
since the zeroes are not easily determined in that case when the lattice is finite, they termed 
the argument as "hand- waving" . Here, the argument is made rigorous by the use of the 
Brascamp-Kunz boundary condition. From here on, therefore. We shall adopt this simpler 
approach in all subsequent considerations. 

We now proceed to determine the density of the zero distribution. Let the number of 
zeroes in the interval [a, a + da] be 2MN g{a)da such that 

/ g{a)da = 1 (8) 
Jo 

Hz-e""). (9) 
Jq g{x)dx where 2MNR{a) gives the 

g{a) = ^R{a). (10) 



/ = - ln(4z) + / da g{a) 
2 Jo 

It is more convenient to consider the function R{a) = 
total number of zeroes in the interval [0, a] such that 



On the circle |2;| = 1 writing z = e" and setting the argument of the logarithm in the 
third line of (^) equal to zero, we find a determined by 

COSa = COSMCOSf , 0<{u,v}<7i. (11) 

Now if ai is a solution, so are —at and vr — a^, hence we have the symmetry 

fi'(") = 5'(-") = 5'(7r - «). (12) 

It is therefore sufficient to consider only < {a,u,v} < 7r/2. 

The constant-a contours of |lT| are constructed in Fig. 1(a) and are seen to be symmetric 
about the lines u,v = ±7r/2 in each of the 4 quadrants. Now from (|^) we see that zeroes 
are distributed uniformly in the {9, cp}-, and hence the {u, f }-plane. It follows that R{a) is 
precisely the area of the region 

cos a > cos M cost), < {a, u, f } < 7r/2, (13) 
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(a) (b) 

Figure 1: Constant-a contours in the u-v plane, (a). The contour ([TTl) for the simple- quartic 
lattice. Straight lines correspond to a = tt/2. (b). The contour (|23|) for the triangular lattice. 
Broken lines correspond to a = 2cos^^(l/3). 



normalized to R{'k/2) = 1/4. This leads to the expression 



1 r 
vr^ Jo 



cos 



cos a 



cosx 



dx. 



(14) 



Using ( pUD and after some reduction, we obtain the following explicit expression for the 
density of zeroes, 

' sin a I 



9W 



R'{a) 



-K{sin a), 



(15) 



where K{k) = Jq^"^ dt{l — k"^ sin^ t)~^/^ is the complete elliptic integral of the first kind. The 
density (|T5|), which possesses an unexpected logarithmic divergence at a = ±7r/2, is plotted 
in Fig. 2(a). For small a, we have g{a) ~ \a\/27T. As pointed out by Fisher it is this 
linear behavior at small a which leads to the logarithmic divergence of the specific heat. 

We can also deduce the density of zeroes on the two Fisher circles (|l]) which we write as 

tanhK ±1 = V2e''^. (16) 

The angles a and 9 are related by, 

•72 Te-'" 



so that the mapping from a to 6' is 1 to 2. This leads to the result 

9(0) 





da 


2 


dO 



(17) 



;i8) 



Let the density of zeroes be g±{6) for the two circles (|T6|). Then, using (^) we find 

K{k), (19) 





1 - V2cos^ 




3 - 2^2 cos ^ 
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where 

* = Si!!i^M^^. (20) 

The density ( [T9|) is plotted as Fig. 2(b). Note that the divergence in the density distribution 
in (p^ ) is removed in (p!9|). The points a = ±7r/2 in g{a) is mapped onto the points 
e = ±ln/2± tt/A) in g±{e). We have g+in/A) = g-iSn/A) = 0, and for small 6 we find 

(21) 

Here, again, the linear behavior of g+{6) at = leads to the logarithmic singularity of the 
specific heat. 

It is also of interest to consider zeroes of the Ising model in the Potts variable x = 
(e^^ — l)/\/2. In the complex x plane it is known ^ that the partition function zeroes are 
on two unit circles centered at x = 1 and x = —\/2. We find the density along the two 
circles to be, respectively, g^{0) and g+{0). 




3 The triangular lattice 



For the triangular Ising model with nearest-neighbor interactions K, the free energy assumes 
the form B |9| 



/ 



C + 



C + 



dO 



>ln 



27r2 Jq 



du / dv In 







z + z + 1 — [cos 6 + cos + cos(^ + 
z + z'^^ +2 — 2 cos u(cos u + cos v) 



(22) 



where C = [\n{Az)]/2, z = (e — l)/2, and we have introduced variables u = {6 + (f)) /2, v = 
{6 — 0)/2. Now the value of the sum of the three cosines in (|2^) lies between —3/2 and 3. It 
then follows from Corollary 1 that in the complex z plane the zeroes lie on the union of the 
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unit circle \z\ = 1 and the line segment [—2, —1/2] of the real axis, a result first obtained by 
Stephenson and Couzens 0. 

The density of the zero distribution can now be computed in the same manner as de- 
scribed in the preceding section. For z on the unit circle we write z = e*". Then a is 
determined by 

cosa = — 1 + cosm(cosm + cosf ), 0<{u,v}<'it, (23) 

and R{a) is the area of the region 

cosa > — 1 + cosm(cosm + cosf ). (24) 

Clearly, we have the symmetry (y'cir(«) = gciriT^ — d) and we need only to consider < a < vr. 

From a consideration of the constant-a contours of ^ shown in Fig. 1(b), we obtain 
after some algebra the result 

^cir(a) = 4^^i^(A:), (25) 
where A{a) = (5 + 4 cos a) and 

= F[A{a)] 

F{x) ^ ^{--l){l + xf. (26) 
lo ' 

Particularly, for small a, we find gdxip) ~ |a|/2-\/37r. 

In a similar fashion we find, on the line segment z G [—2, —1/2], we write z = —e^ and 
obtain 

g^^^{\) = A^^^^K{k-^), -ln2<A<ln2, (27) 
where -B(A) = [5 — 4 cosh A]^/^ and 

= F[B{\)\. (28) 

In contrast to the case of the simple-quartic lattice, the density of zeroes is everywhere finite. 
Specifically, we have 5'cir(7r) = 5'ime(0) = 0, and 5'iine(±ln2) = v^/27r. The densities (pSf ) and 
(^) are plotted in Fig. 3. 

Matveev and Shrock have discussed zeroes of the triangular Ising model in the 
complex u = e~'^^ plane, for which the zeroes are distributed on the union of the circle 

7/ = l(2e^" - 1), -7r<a<7r, (29) 

and the line segment 

-oo<M<-^. (30) 
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Figure 3: Density of partition function zeroes for the triangular lattice, (a), gdrici) given by 
(H). (b). g\ine{z) given by (|^. 



Using our results we find the respective densities 

/ N I sin 



97r2 



where C{a) = 3(5 - 4 cos 0)-^/^ = F[C{a)], and 

(1 + m)(1 -3m) 



9\me{u) 



AT^'^u'^il-uYk^ Diu 



Kik~'), 



(31) 



(32) 



where D{u) = y (1 + 3u)/u{l — u) and /c^ = F[D(m)]. At the end point we have 5'iine(— 1/3) = 

The density of zeroes assumes a simpler form if we use Corollary 2 to map all zeroes onto 
a unit circle in the complex w plane, where w is root of the quadratic equation 



(33) 



and z = {e^^ — l)/2. For w on the unit circle, we write w = e*" and find in analogous to 
(|T3|) that R{a) is the area of the region 

1 



cos a > 



9 



1 cos m(cOS U + COS v) — 7 



(34) 



Using the contours shown in Fig. 1(b), we obtain 



R{a) 



TC"^ Jo 



00 

cos' 

1 



9 cos a + 7 
8 cosd) 



cos ( 



2 TT" J<j,o 

where ao = 2 cos~-'^(l/3) and 

(j)r) = cos~^ 



IT — COS 



COS 



9 cos a + 7 
8 cosd) 



COSi 



a E [0, ao] 
a e [ao,7r]. 



(35) 



3 


a 


1 




cos — 




2 


2 


~ 2 




i[3„„ 


a 



a 1 
tor cos — < - . 



(36) 
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Note that we have R{ao) = 3/8, R{n) = 1/2. 
Finally, using (p!OD, we obtain 



9 sin a r<t>o 



iTC^ Jo 



9 sin a r<t>^ 



'cos^ (p — cos^ 0o)(A2 — COS^ I 



^ V (COS^ — COS^ 0o) (COS^ 01 — cos^ 



a G [0, ao] 
=, a G [ao,vr] 



(37) 



where A = [1 + 3 cos(q;/2)]/2. After some manipulation and making use of integral identities 
l1|) and ([A2D derived in the Appendix, we obtain 



3^3 



87r2 
3^3 



87r2 



sm a 
sin a 



sec —K{k), a G [0, ao] 



sec — ^), aG[ao,vr] 



where 



a 



/c = — (sec — - 1)( 1 + 3cos — 
lo 2 



a\3 



(38) 



(39) 



Note that g{a) diverges logarithmically at a = ±ao. 



4 Simple- quart ic Ising model in a field in/2 

The two-dimensional Ising model can be solved when there is an external magnetic field 
z7r/2. The solution for the simple-quartic lattice was first given by Lee and Yang |Ty] and a 



rigorous derivation of which was given later by McCoy and Wu In 1988 Lin and Wu 
||r^ gave a general prescription for writing down the solution of the Ising model in a field 
in/ 2 by transcribing the solution in a zero field. The most general known solution of the 
Ising model in a field 27r/2 is a model with a generalized checkerboard type interactions [0. 
We consider in this section the case of the simple-quartic lattice. 

For the simple-quartic lattice Lee and Yang [TO] gave the free energy in a field in/ 2 as 



/ = i- + C + ^-^ r de r d01n[z + 2-i + 2 - 4cos^cos0], (40) 

where C = (Insinh 2iir)/2, ^ = e"^^ . Setting the argument of the logarithm in (^Op equal 
to zero we have —6 < z + < 2 and hence from Corollary 2 we see that in the complex 
z plane zeroes of the partition function lie on the unit circle |z| = 1 and the line segment 
—3 — 2\/2 < z < —3 + 2\/2 of the real axis. 

On the unit circle |z| = 1 we write z = e*" and find the density 

9cUa) = ^-^K{k), (41) 



2tt^ 



where 



A;2 = (3 + cos a) (1 - cos a) /4. (42) 



8 



On the line segment, we write z = -e^ with -2 ln(l + v^) < A < 2 ln(l + ^2), we find the 
density 

^Hnc(A) = ^lJ^i^(fc), (43) 



where 



fc2 = (3-coshA)(l + coshA)/4. (44) 

At the end points we have 5'iine(±2 ln(l + v^)) = 1/V2n. The density functions (^Tf ) and 
(H) are plotted in Fig. 4. 
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Figure 4: Density of partition function zeroes for the simple-quartic lattice Ising model in a 
field i7i/2. (a). 5'cir(a) given by (M). (b). 5'iine(^) given by (H). 



5 Triangular Ising model in a field 27r/2 



The solution for the triangular model in a field in/ 2 was first obtained in |]T2| by applying 
a transformation in conjunction with the solution of a staggered 8-vertex model. Here, for 
completeness, we present an alternate and more direct derivation. 





(a) 



(b) 



Figure 5: (a). The triangular lattice, (b). A unit cell. 
Consider a triangular Ising lattice of N sites whose sites are arranged as shown in Fi^ 
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5(a). After making use of the identity e 



ia, the partition function assumes the form 



E n 

(Ti=±l nn 



J<CJiCJj 



n 



(45) 



where the first product is over all nearest neighbors, and the second product over all sites. 
Now it is known that the triangular Ising model can be mapped into an 8-vertex model on 
the dual of the square lattice also of sites. However, in order to properly treat the 
factor Hj in (^5]), we need to divide the N "cells" of the lattice, where a cell is shown in 
Fig. 5(b), into two sublattices, A and -B, and associate two cr^'s to each cell belonging to 
one sublattice, say, B. This permits us to rewrite (|45|) as 



■N 



E n ^'^tg(^l'^2,0-3,0-4) 
"i=±l cells 



(46) 



where Wstg{cri, (72, cs, 0-4) is a staggered Boltzmann weight given by 



W^stg(0"l, 0-2, 0-3, 0"4) 



((Tl (72 +(T2 0-3 +CT3 (Tl ) 



for A 



Icri0-2je 



^(cti 0-2+0-20-3+0-30-1) ^ 



The 8-vertex weights are 

{ui,---,ujs} = {e 



3K 



3K 



-K 



-K 



-K ^iK^ 



-K 



-K 



-K 



-K 



-K 



-K 



}■ 



(47) 



(4J 



Furthermore, from the mapping convention of Fig. 1 of |]I5[, we see that the mapping between 
the spin and 8-vertex configurations is 2 to 1. This leads to 



2z^Z;v({^},M), 



(49) 



which is an exact equivalence between and the partition function Zn{{uj}, {uj'}) of the 
staggered 8-vertex model. 

Now the weights ( ^8|) satisfy the free-fermion condition [jl4| for which Z]^{{uj}, {oj'}) has 
already been evaluated |jl5|. Using Eq. (19) of [0 and after some reduction, one obtains 
the following expression for the per-site free energy. 



/ 



c + 



de 



>ln 



1 + e^^f + 4cos0(cos6' + cos( 



2 ■ " ■ 47r2 

where C = [ln(2 sinh 2ii')]/2. As a result, the partition function zeroes are located at 



[l + e 



iK\2 



—4 COS </)(cos 9 + COS ( 



o<{e,0}< 



TT. 



It is therefore convenient to consider the z = e plane. Since 

-8< (l + e^^)2<l. 



(50) 



(51) 



f52) 



using the Lemma we find that the zeroes are on the union of the segment —2 < z < of the 
real axis and the line segment z = —1 + iy, —2^2 <y< 2\/2. The density of zeroes can be 
similarly determined. On the segment z G [—2, 0] of the real axis, we find 



'l + z) 



27T^k^ E(z 



K{k 



-1^ 



(53) 
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where E{z) = \l-z{2 + z) and k"^ = F[E{z)]. Particularly, we have g{0) = g{-2) = l/V^ir 
and — 1) = 0. On the line segment z = —1 + iy, we find 



9iy) 



y 



2n\Hiy) 



Kik), 



(54) 



where H(y) = y/1 + y^ and k'^ = F[H(y)]. Particularly, we have g(0) = and g{±2^/2) 
l/-\/67r. These results are plotted in Fig. 6. 
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Figure 6: Density of partition function zeroes for the triangular Ising model in a field z7r/2. 



(a), giz) given by (|3]). (b). g{y) given by 

we remark that in the complex x = e~^^ plane considered in |16], the segment —2 < -2 < 
of the real axis maps onto — oo < x < —1/2 while the line segment z = —I + iy, —2^/2 < 



y < /y/, IS mapped onto the circular arc x = 
The density of zeroes on the arc is found to be 

(1 + cos 0)2 

5'arc(6') 



U-1 + e'^), 00 = tan-i(4v^/7) < \e\ < tt. 



2tt\ i{e)\sm'e\ 



Kik), 



(55) 



where 1(9) = [2(1 + cos0)]^/V| sin 6\ and k'^ = F[B{6)]. The densities at the end points of 
the arc are gard^do) = 3V3/2^/2tc. 



6 The honeycomb and kagome lattices 

The partition function of an Ising model on a planar lattice with interactions K is propor- 



tional to the partition function on the dual lattice with interactions K* jl^, where K and 
K* are related by 

= tanhir. (56) 

Consequently, their partition function zeroes coincide when expressed in terms of appropriate 
variables. Now the honeycomb and triangular lattices are mutually dual, it follows that for 
the honeycomb lattice with interactions K, in the complex 

z = he^^* - 1) = (cosh2ir - 1)-^ (57) 
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plane, zeroes of the partition function coincides with those of the triangular lattice partition 
function (j2|). 

For the honeycomb Ising model in an external field z7r/2, the free energy can be obtained 
from that in a zero field via a simple transformation [0, . Writing the partition function 
in the form of (^Sf) and replacing the product YiiCi by Hjcf; it is clear that, besides the 



factor , the partition function is the same as that in a zero field with the replacement 

or, equivalently, e^^ — > —e^^. It follows that in the complex 

z = {- cos\i2K - l)-^ (59) 

plane, the zeroes coincide with those of the triangular lattice partition function (p^). 

The Ising model on the kagome lattice with interactions K can be mapped to that on an 
honeycomb lattice with interactions J, by applying a star-triangle transformation followed 



by a spin decimation. The procedure, which is standard |[T^ and will not be repeated here, 
leads to the relation 

g2J ^ ^g4i^ ^ ^y2. (60) 

As a result, we conclude that, in the complex 

z = (cosh2J- 1)-^ = 2(1 -tanh2K)/tanh2 2Js: (61) 

plane, zeroes of the kagome partition function coincides with those of the triangular lattice 
partition function (0). The evaluation of the kagome partition function in an external field 
in/ 2 remains unresolved, however. 
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APPENDIX Two integration identities 

In this Appendix we derive the integration identities 



J, . r , =-^kL ^] (ai) 

•^0 7(1 -a2 smH){b^ + a^ sin^ t) Va^T^ \ ^ a'' + by 



b dx ^ 1 1 b^-a^ \ 

^(1 - x2)(x2 - a2)(62 - a;2) ~ 6^1 - a2 VfeVl-aV' 

which do not appear to have previously been given. 



To obtain (|Al|), we expand the integrand using the binomial expansion 



;i - X)-" = f (A3) 

fc=0 
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where {ct)k = Q;(q; + 1) ■ ■ ■ (a + A; — 1) = r{a + k)/r{a), and carry out the integration term 
by term using the formula 

- / sin2™tdt = ili^. (A4) 
vr JO ml 



This yields 



IT 



^ ^■IHi'a- _L i-^l V 1)2/ 



26,t^o j!A;!(j + A;)! 



TT 



1 1 1 



(A5) 



where the hypergeometric function of two variables is (cf. 9.180.1 of [19|) 



(A6) 



This leads to the integration formula (|A1|) after making use of the identity (cf. 9.182.1 of 

F,{a; (3, (3 + (3'; x, = (1 - y^-F (^a, f3; (3 + (3'; , (A7) 

where F is the hypergeometric function (cf. 9.100 of |]19 



Fia,P;r.z) = t^f^z^ 



(AJ 



and the identity (cf. 8.113.1 of p|) 



(A9) 



The integral ( |A^ ) is obtained by introducing the change of variable x"^ = (62 — a2) sin^ t + 
a^, which yields 



1 



1 — a2 Jo 



tt/2 



dt 



1 - c2 sin2 1) [a2/(l - a2) + c2 sin^ t) 



(AlO) 



where = [b"^ — a^)/(l — a^). The integral I2 is now of the form of Ji and ( [A^ ) is obtained 
after applying ( [Al| ). 
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